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1.  INTRODUCTION  AND  SUMMARY 

This  report  describes  the  fluctuations  In  intensity  that  occur  when 
a radio  signal  traverses  striated  ionization.  The  ionization  is  modeled  as 
a collection  of  cylinders  of  electron  density  randomly  positioned  in  space. 
This  work  is  a continuation  of  a previous  analysis^.  The  major  in^rovements 
on  that  analysis  are: 


1.  The  use  of  fast  Fourier  transforms  to  calculate  the  signal 
a distance  d,  beyond  the  ionization  instead  of  the  method 
of  rays  and  caustics  used  previously. 

2.  The  analysis  is  carried  out  for  two  extreme  profiles  of 
striation  electron  density;  the  Gaussian  profile  previously 
studied  and  a rod  profile  representing  an  extreme  steepening 
of  the  striation  edge. 

3.  Results  of  systens  significance  are  presented  graphically 
as  contours  on  a Xd  versus  plot  where  X is  the  signal 
wavelength  and  is  the  root-mean-square  signal  phase 
variation  directly  produced  by  the  ionization. 

The  last  i^rovement  allows  the  difference  between  the  extreme  profiles'  effects 


to  be  easily  se«i  as  a function  of  the  parai^ters  Xd  and  A* 


rms 


.1 


The  set  of  striations  used  is  the  s»ne  as  that  used  previously  , a set 
of  210  striations  set  randomly  in  an  area  64  km  wide  by  a thickness  T.  The 
thickness  does  not  enter  the  calculations  since  use  the  thin-screen 
approximation.  This  approximation,  v^ose  validity  criteria  were  previously 
derived^  for  Gaussian  profiles,  is  adequate  to  d^nstrate  the  propagation 
effects  for  the  striations  considered  in  this  report.  We  do  not  foresee  any 
qualitative  changes  by  its  rswoval  and  intend  to  demonstrate  this  in  future 
works. 


We  find  the  rod  and  Gaussian  profiles  to  produce  similar  propagation 
2 

effects  for  Ad  > 3 km  where  the  dependence  on  Ad  is  weakening.  In  the 


high  Ad  regime,  amplitude  scintillation  onsets  at 


0.3  radians. 


As  Ad  becomes  smaller  the  rod  and  Gaussian  profile  effects  diverge  with  onset 
for  Gaussian  effects  occurring  at 


/.A 

’^nns 


whereas 


for  the  onset  of  rod  effects. 

-3 

Thus  the  profile  shape  is  relatively  unimportant  at  UHF  (A  10  km), 

2 4 

but  becanes  very  important  at  k Band  (A  -v  10  to).  We  ass»m«  10  < d < 10  to 

3 

for  typical  situations.  For  d » 10  and  K Band  the  Gaussian  profile  effects 
onset 


at  = 100  radians 

rms 

while  the  rod  profile  effects  onset 


at  * 6 radians, 

ms 

The  rms  phase  is  given  by 


"^%s  * 2-818.10“^A 

“■3 

vdiere  n^  is  the  striation  on-axis  electron  density  (cm  ) and  Q is  the  normalized 


4 


integrated  electron  density.  For  the  particular  set  of  striations  used, 

^nns  “ km  for  the  Gaussian  profile  and  = 2.55  km  for  the  rod  profile. 

This  indicates  Av^  to  be  relatively  insensitive  to  profile  shape  for  a 

given  set  of  striations  and  values  of  Oa  and  X.  However,  the  value  of  A# 

required  for  the  onset  of  degradation  is  highly  sensitive  to  profile  when 
2 

Xd  < 3 km  . The  set  of  striations  we  have  used  have  sizes  ranging  from  323m 

to  I61te  with  a maximum  mmrt)er  at  775m.  One  can  scale  our  results  to 

*■2 

different  sizes  by  changing  each  size  a to  5a  and  scaling  Xz  as  5 • The  size 
distribution  used  is  a discrete  approxi#at1on  to  the  continuous  Chesnut 
distribution  which  gives  the  analytic  formula  for  Gaussian  profiles 


3 


where  3q  is  a characteristic  size  of  the  distribution  {'''  0.76  km  for  the  set 
used),  M is  the  total  nwnber  of  striations  imbedded  in  an  area  L x T.  Since  T 
does  not  matter,  if  we  scale  a^  and  L to  and  6L,  scales  as  5. 

This  study  has  contrasted  the  propagation  effects  of  smooth  versus 
abrupt  profiles  for  a set  of  striations.  The  abrupt  profiles  produce  more  severe 
propagation  effects  than  snwoth  profiles  at  high  frequencies.  Future  w>rk  will 
address  this  question  for  other  sets  of  striation  sizes.  We  do  not  foresee 
other  sets  of  striations  producing  effects  more  severe  than  were  found  here 
for  the  abrupt  profiles.  We  illustrate  this  statement  with  an  exar^le. 

Suppose  we  are  at  the  threshold  of  rod  effects  at  Xd  = 0.1  and  ^^rms  ' 

If  we  change  the  sizes  of  all  striations  by  dividing  then  and  the  width  L by 
y/W we  effectively  have  changed  Xd  to  1 and  to  1 which  keeps  us  at  the 
threshold  of  rod-effects  (Gaussian  effects  get  worse).  If  we  continued  to 


5 


divide  the  sizes  we  would  reach  < 0.3  and  propagation  effects  would 

cease.  We  believe  the  results  presented  in  this  report  for  rod  profiles 
constitute  a reasonable  upper  bound  for  estimates  of  system  degradation  e 
In  terms  of  and  Xd. 


2.  ANALYTIC  METHODS 


2.1  ANALYSIS 

We  use  the  thin-screen  approximation  to  describe  the  signal  which 
results  after  a plane  wave  propagates  through  a region  of  striated  ionization. 
The  derivation  and  validity  criteria  of  this  approxiTOtion  are  in  reference  1. 
The  ionization  is  represented  by  a collection  of  striations  each  of  which 
has  a profile  of  electron  density 


OgCx)  = n.{x-x^,a.) 


where  is  the  two  di^nsional  position  of  the  i‘th  striation  and 

is  its  characteristic  size.  The  electron  density  is  uniform  in  the  third 
dimnsion. 

For  a plane  wave  propagating  In  the  z direction,  the  thtn-screen 
approximation  describes  the  effect  of  each  striation  on  tte  signal  to  be  a 
phase  change 


a^)  ® " *i» 


where  kQ  is  the  wave  nwober  of  the  signal  = 2»/X  where  x is  the  wavelength 

1 -13 

and  n^  * •'o  critical  electron  density  where  = 2.818.10  m*  the 

electron  classical  radius.  The  total  effect  on  the  signal  is  then  a net  phase 

change  doe  to  all  the  striations  of 


#Cx)  = E#|Cx  - X.,  a|)  . 


He  consider  two  extremes  of  electron  density  profile : 


2,.  2 


The  Sausslan,  n^S  = ~ 


tAich  produces  the  phase 


2,.  2 


A . 


9 = 


‘{x  “ X^}  /d* 


the  rod.  n 


* — ah  1*^  ( 

i ■ "o  •**  " ’‘il  ^ ®i 


0 !x  - x.f  > 8| 


idiich  produces  the  phase 


r / C*  * 

■ “•n/' 


'I 


where  C . • V*0  • 

c 


The  actual  profile  of  a striation  stould  fall  between  tNise  tuo 


extremes. 


The  total  phase  has  the  Fourier  trMsform 


F(v) 


e'^ 


i2svx 


‘f(x) 


Ee  ' 


»Aere  F.®(v) 


F|Cv) 

J 2,  2 

2 -s  V a 


Cta^^  e ’ and  F.’’(v)  ^ Jj(2sva|), 


w^re  J,  is  a Bessel  fimtiMi. 


Notice  that  F^-^(O)  = F^^(O)  = Cna^^  because  each  profile  has  the  same  electron 
content,  /ix/diZ  n^(x  - x^,  a^.)  = unga,^. 


For  a given  set  of  striations,  the  phase  has  a power  spectrum 


P(v)  = lF(v)l^  = Lf.^(v)  + F^.(v)F.(v)  cos[2ttv(x^.  - x.)]  . 


The  first  term  is  independent  of  the  positions  of  the  striations.  The  second 


term  would  be  larger  than  the  first  by  about  a factor  of  the  total  number 


of  the  striations  if  the  cosine  terms  all  added  coherently.  The  discussion  of 


the  propagation  effects  of  a collection  of  striations  should  ideally  not  be 


dependent  on  their  particular  positioning  unless  it  can  be  shown  that 


characteristics  of  their  positioning  can  be  predicted  by  the  phenomenology. 


Since  no  positioning  characteristics  have  yet  come  from  pheonomenology,  the 


second  term  can  be  eliminated  by  considering  the  striations  to  be  randomly 


placed  and  averaging  over  an  ensemble  of  random  placements.  If  the  total  area 


occupied  by  the  striations  is  small  compared  to  the  area  they  are  placed  in, 


one  can  approximate*^  the  probability  of  any  two  striations'  placement  as 


independent  and  equally  likely  within  the  area.  The  position,  x^,  is  then 


equally  likely  on  a line,  L,  which  produces  an  average  power  spectrum  of 


^ + 2S^  L F.(v)F.(v) 

-*  J 4^4  * J 


Where  . 


By  evaluating  the  spectrum  at  values  of 


-p  n = 1,2,3. 


where  S = 0,  we  are  correctly  ignoring  the  effects  of  the  finite  interval, 

L,  on  the  Fourier  transform.  We  are  expanding  $(x)  in  a Fourier  series. 

The  propagation  calculations  will  be  performed  for  one  set  of  random 
positions  of  the  striations.  This  is  the  situation  an  actual  radio  signal 
encounters.  We  illustrate  in  Figures  1 and  2 the  power  spectrum  for  the 
one  set  and  the  average  power  spectrum  for  the  Gaussian  and  rod  profiles. 

The  scale  is  normalized  (C  = 1).  We  use  a set  of  positions  and  radii  furnished 
by  Dr.  W.  Chesnut  (SRI).  This  is  his  "Type  B"  distribution  consisting  of 
210  striations  whose  positions  are  randomly  distributed  over  a line  L = 64  km. 

This  distribution  was  also  used  in  reference  1. 

It  is  obvious  that  the  large  differences  between  the  specific  set  of 
positions  and  the  average  at  low-v  is  a chance  alignment  in  phases  for  the 
positions  chosen.  It  appears  for  both  Gaussian  and  rod  spectra. 

The  rod  spectra  have  a higher  frequency  content  of  course.  The  spectra 
appear  to  change  from  exponential  to  power  law  at  v = 1.  Since  there  are  no 
striations  in  the  set  with  a^  < 0.323  km,  all  Bessel  functions  would  have 
arguments  greater  than  2 for  v = 1.  They  then  decrease  with  an  envelope  of 

giving  the  power  spectrum  a v'^  fall  off.  The  Gaussian  spectrum  is 
exponential  until  v = 1,5  where  the  Gaussian  fall  off  should  take  over 

The  propagation  calculations  we  will  perform  with  the  specific  set  of 
striations  depend  on  the  Fourier  transform  of  We  cannot  yet  analytically 

determine  the  effect  of  the  large  low-v  variation  in  the  power  spectrum  of 
♦(x)  for  the  specific  set  on  our  calculated  effects.  Our  intuition  says  the 
low-v  portion  of  the  spectrum  is  unimportant  so  long  as  it  is  finite.  We  also 
feel  that  all  statistical  properties  of  the  signal  of  interest  will  be 
unaffected  by  the  noise-like  difference  between  the  specific  and  average  spectrum. 
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Consider  the  following  statistical  functions: 


III  5 

f = -J_  f:  FifjeiS’-VL  . 


n=l  ^ >- 


That  iSj  the  Fourier  transform  of  A(j>  is  the  same  as  the  Fourier  transform  of 
4>  with  the  v = 0 term  removed.  The  auto  correlation  function 


y*dx  A<{)(x)A(jj{x  + Ax) 


p(0)  fei2TrnAx/L  ^ g.i2nnAx/L  3 

n=l 


which  follows  by  substitution  and  the  properties 


0 n f 

1 n = -n'  , 


F(^)  = F*{f)  . 


We  have  the  well  known  result  that  the  power  spectrum  is  the  Fourier 
transform  of  the  auto  ccrrelation  function.  In  this  case  we  are  referring  to 
power  per  unit  length  since  we  are  interested  in  A(|?  over  the  interval  (0,L)  only. 
Over  this  interval  the  auto  correlation  function  is  meaningful  out  to  |ax|  = L/2 
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since  beyond  this  value  the  periodicity  of  A#  causes  it  to  mirror  its  behavior 
back  to  its  value  at  Ax  = 0.  That  is,  /V^Tx^AoIx+l-Ax ) = aT(x)Ac (x+AxT. 

Indeed,  for  cases  of  interest  to  be  treated  later,  = -^Z^PCr)  » 1 

n ^ 

This  causes  most  of  the  interest  in  the  auto  correlation  function  to  be  in 
its  behavior  at  small  ax.  For  example,  one  item  of  interest  is  that  value 
of  AX  where 


— 7 

A$A(|)(ax)  - hfy  = -1  or 


L2EP(t)  n - cos(2Tt^^^)]  = 1 or 


4EP(f^sin2(-2^)  = 1 
in'- 


Since  this  will  occur  at  a value  of  Ax  « L,  the  low  v (small  n) 
difference  between  the  specific  and  average  spectrum  should  not  impact  the 
result.  Similarly  for  ttAx  « L,  the  sin  weighting  will  vary  slowly  over  a 
range  of  n values  which  smoothes  out  the  differences  between  the  specific  arc 
average  spectrum. 

The  discrete  distribution  of  that  we  use  approximates  the  connnuous 
Chesnut  distribution. 


m(a)da 


8M  ®o 

— — r-  e ' 0 'da  , 
a® 


the  number  of  striations  of  size  a in  an  interval  da.  The  total  number  of 
striations,  M = 210. 


Analytically, 


2 2 2 

PTvT  = _/in(a)da  |c"a^  e'^  ^ ® j 


_ 4 f.2„2,  4„ 

- ^ C . M e 


for  the  Gaussian  profile. 

Figure  1 shows  the  exponential  behavior  quite  well  except  near  v = 0 
where  the  cutoff  of  large  striations  in  the  discrete  distribution  returns 
the  vanishing  slope  at  the  origin.  From  Figure  1,  a measured  slope  of  7.83 
produces  the  value  a^  = 0.881  km. 

We  are  interested  in  relating  the  propagation  effects  to  a character- 
istic size  of  the  distribution  used  and  to  the  parameter 


f cV(") 

L n n 


yfB  3/Mv 


for  the  continuous  distribution.  We  have  measured  this  quantity  directly 
from  the  discrete  distribution  and  find 


^ = 2.071 1C 


Substitution  of  this  value  into  the  abov2  equation  produces  = 0.761  km. 
Thus  there  is  a 10%  discrepancy  in  determining  from  specific  and  average 
discrete  distributions  using  formulas  for  a <ontinuous  distribution.  Since 
the  geometric  structure  of  the  striations  we  use  is  fixed,  we  will  relate  the 
propagation  effects  to  A?rms  variations  of  the  parameter  C. 
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We  now  discuss  the  signal  at  a distance  d beyond  the  striated  region. 


If  the  signal  Is  defined  as 


E = U{x.z)e 


l(k^z-wt) 


and 


U(x,0)  = at  the  exit  plane  of  t*i?  '.tri.":.ions,  the  signal  at 


some 


point  (xQ,d)  beyond  the  exit  plane  is  cbUir  d f 'he  free-space  parabolic 
equation 


3U  . 1 3 U _ n 

5?  - » • 


Fourier  analysis  of  the  parabolic  equation  shows  that  for 


U(x,0)  = j;u(k)e*'‘’‘ 
k 


k^d 


^ iCkXo  - Jjr  3 

u(x  ,d)  = £u(k)e  ® 

° k 


This  method  has  had  widespread  us^  by  workers  in  striation  propagation  and 
is  efficient  for  production  runs  of  U(xQ,d)  because  of  fast  Fourier  transform 
techniques.  This  technique  decon^ses  the  signal  at  z * 0 into  plane  waves, 
propagates  each  plane  wave  to  the  receiver  at  (x^jd)  and  sws  the  waves  to 
obtain  U(xQ,d).  Since  U{x,0)  was  defined  only  for  0 < x < L,  the  use  of  a 
sum  {Fourier  series)  rather  than  an  integral  specifies  U(x,0)  as  periodic  in  L. 
We  can  accept  this  as  an  unin^ortant  consequence  so  long  as  we  do  not  evaluate 
the  signal  at  values  of  d so  far  from  the  screen  that  the  signal  is  influenced 
by  values  of  U(x,0)  over  an  interval  greater  than  L. 
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The  Fresnel  distance  for  an  interval  L is 


= kgLV2  . 


3-1  7 

The  smallest  wave  number  of  interest  is  -vS-lO  km  so  we  need  d < d^  w 10  km 
for  L = 64  km,  which  is  well  satisfied  for  cases  of  interest.  A further 
requirement  when  large  phase  variations  occur  is  L/2d  > 6^  where  8^^  is  the 
largest  angle  of  a significantly  contributing  plane  wave. 

The  parabolic  equation  is  a small  angle  approximation  to  the  wave 
equation.  For  each  plane  wave,  the  parameter  k is  k^^sinej^  ^o®k  ' 
amplitude  of  the  plane  wave  is  the  Fourier  integral 


^g1[<l'(x)-kx]^j^ 


As  in  a previous  analysis  , we  may  define  a local  phase-front  angle 


0(x)  * IP  H 


and  the  integral  is 


U(k)  = -T“  i ® ' 


If  instead  of  fast  Fourier  techniques,  we  evaluated  the  integral  by  stationary 
phase^,  the  points  of  stationary  phase  wuld  be  wherever 


eCXg) 


and 


fl(ic)  = E 


itO(,(Xs)-kX5] 


with  a criterion  for  the  validity  of  this  method  of  ]0"(x  )j^  « k J5'(x 

Although  this  method  of  evaluation  of  the  integral  becomes  inaccurate 
where  8"(x)  is  small,  it  is  clear  that  there  are  negligible  contributions  to 
the  integral  for  6.  > since  there  would  then  be  no  stationary  points. 

K fiiaX 

For  a Gaussian  profile 

-2f -r  -(x-x.)^/a.^ 

6i(x)  = (x-x.l  e ’ ' 

’ Vi  ' 


with  maximum  value 


e K 


Iq  n. 


and  depending  on  the  number  of  striations  along  the  path  there  is  some  finite 
limit  of  9j^  for  which  the  requirement 


can  be  tested. 


For  the  case  of  a rod  profile 


0^(x) 


(x-X.) 

"”(x-x^ 

' “T2“ 


which  approaches  infinity  as  !(x-x^-)i  approaches  a^  and  there  is  no  limit 
on  0(x).  Indeed  there  are  two  infinities  per  striation  in  6(x). 


However 


■3^{x) 


C(x-x.)2/a.2  + 1 


(x-x-K 


also  becomes  infinite  so  the  amplitude  of  the  plane  wave  contribution  goes 
to  zero. 

The  rod  profile  is  an  extreme  idealization  of  the  sharpness  of  the 
edge  of  a striation.  The  fact  that  we  will  evaluate  U(k)  by  a fast  Fourier 
transform  rather  than  analytically,  means  we  are  automatically  smoothing  the 
edges  of  the  striation  since  the  largest  value  of  6 obtainable  is 


;.{x  = x.+a.-Ax)  - <{).(x  = x^+a^) 


n /2*a7 

So  ^ 7-  \t  where  Ax  is  the  step  size  = L/N  where  N = 2"^  is  the 
\rax  \ V 

number  of  Fou»*iGr  components  (plane  waves)  to  be  determined. 

We  have  carried  out  the  analysis  for  m = 14  and  15  to  see  the  impact 
on  prcoagation  effects  of  changing  the  step  size  and  hence  the  maximum  plane 
wave  angle,  "he  differences  were  negligible,  indicating  that  for  the  cases  we 
■■.ave  treated,  the  largest  angle  of  a significantly  contributing  plane  wave  was 
less  t"an  the  largest  angle  obtainable  by  using  m * 15.  We  do  not  know  what 
this  angle  is  for  the  rods.  We  estimate  the  angle  for  the  Gaussian  profile 
using  previous  analysis  which  shows 


If  we  choose  the  restriction 


k„a„L 

o 0 


AA 

■^rms 


we  should  be  confident  that  the  finiteness  of  the  interval,  L,  will  not 
impact  the  calculations.  It  turns  out  that  the  results  of  importance  for  this 
study  occur  well  before  this  limit  is  reached. 

The  parameters  Xd  and  (or  C)  uniquely  determine  the  propagation 
effects  for  a given  striation  structure.  This  can  be  seen  from  an  alternate 
expression  for  the  signal  previously  obtained^ 


U(*o,d)  * 


n~  f i[^^{x)  + ir{x-x  )^/Xd] 


For  a fixed  structure,  A$(x)  and  are  both  linear  in  C and  therefore 
scale  together.  The  only  other  parameter  in  the  integral  is  Xd.  We  may 
therefore  do  all  calculations  at  one  frequency  varying  C and  d,  present  the 

results  in  terms  of  Xd  and  and  obtain  results  at  other  frequencies  by 
scaling. 


2.2 


IMPLEMENTATIW 


We  choose  a range  of  Ad  of  interest  by  considering  the  frequency 
range  of  interest  to  be  300  HHz  to  30  GHz  and  the  distance  range  of  interest 
to  be  10^  to  10^  kffl.  This  limits 

10"2  < Ad  < 10  km^  . 

We  choose  ^/O.l  < ^$^,5  £ 100  radians.  Effects  are  negligible  at  and  below 
the  lower  limit.  The  higher  limit  is  chosen  so  as  to  insure  adequate  sampling 
of  by  the  fast  Fourier  transform.  These  limits  are  sufficient  to 

deronstrate  results.  The  transparency  of  the  mediira  is  assured  for  high 
altitude  plasma  by  the  restriction  £ 100  radians.  For  a temperature 
around  1000®K,  the  ion-electron  moronttm  transfer  collision  frequency 

Vg^{s‘b  **  hg(c«"^)  . 

A wave  traveling  a thickness  2 of  the  plasma  then  suffers 
2 2 

Og  A d8  of  attenuation  where  A is  tf«  wave  length  in  km. 
of  the  striations  the  average  attenuation  is 

I . /dx 

0 0 ’ 


/ d2  1.28-10'^^ 

O'' 

For  the  64  km  width 


Now 


dz  e 


/dx  y 

0 CT 


for  Gaussian  profiles  and 
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r < 


\"a'  " ’il  Z 2 

Jdx  Jdz  I 0 , r > aj  " ="o 
0 0 ’’ 


for  rods. 


Therefore,  the  »-ods  produce  twice  the  average  attenuation  that 
Gaussian  profiles  produce  for  the  saae  electron  content.  For  Gaussian  profiles 
we  have 


A = 10"^^ 

0 *j  1 


For  the  particular  set  of  210  striations  we  use 


= 112.6: 


_ _ii  22  -% 

and  A = 3.54-10  >.  n„  dB  where  X is  in  ka  and  n„  is  ca  . Similarly,  we 

0 0 


= 2.818-10"-^Xn  0 ^ 

rms  o ras 


where  Q{ 


<x)  = 


n^Cx.2) 


is  the  normalised  intearated  electron  der.sitv  in  ktn. 


For  the  particular  set  used,  = 2.071  ta  so  , 

^rm  >Tss 


for  5 100,  A < 10“"  dB.  For  the  rods, 


5,836 -10  "'n. 


y , bs  So 


Ac„  = 7.1S6-10"'Xn  but  A = 7.08'10'^^X^  = 1.37-i0‘®  Ac^  absorpt 
rws  0 0 

is  negligible  for  .t*  ^ < 100.  Our  procedure  is  to  obtain  Qfx)  directly 
the  210  striations  in  steps  of  64  = 1.95te  and  store  it.  For  eacr 
value  of  desired  m miltiply  QC*)  by  to  obtain  :(x)  and 


fast  Fourier  transform  to  obtain  U(k).  We  use  the  Berkeley  computer 

library  routine  CFFT.  Their  write-up  describes  the  inverse  transform  as 


U(x,0)  = 


N-1  . 

S U(j)e 
j=0 


where 


k - M. 
3 L 


This  reproduces  the  function  U(x,0)  exactly  at  the  sampled  points  Xp  = 
p = 0,1, 2... N-1.  Unfortunately,  it  Is  not  the  appropriate  inverse  transform. 
The  proper  one  for  best  fit  to  the  function  at  all  x is 


ik,x  -ik,x 

U(x,0)  = £ U(j')e  £ U(N-j)e  ^ 

j=0  j=l 


which  also  reproduces  the  function  exactly. 

This  is  the  proper  truncation  of  the  infinite  Fourier  series  which  has 
both  positive  and  negative  wave  numbers.  The  one  described  in  the  Berkeley 
write-up  replaced  the  negative  wave  numbers  by  wave  numbers  one  period  higher 
which  introduces  high  frequency  (wave  number)  content  the  original  function 
did  not  have. 

Having  the  U(j)  for  the  specified  we  t! . n modify  U(j)  for  each 

of  a set  of  values  of  Ad  to  describe  the  wave  a distance  d away  as 


-S,.,  -ik/dX/4, 


U(j,d)  = U(j)  e ^ 


j = 0,1,.. N/2  ; 


^ A -ik.'^Xd/4Tr 

U(N-j,d)  = U(N-j)  e j = l,2,..N/2-l  . 


and  call  the  inverse  CFFT  to  obtain  U(x,d).  The  routine  performs  the  inverse 
transform  properly. 

For  the  rod  profiles  we  redid  the  calculations  using  both  2^^  and 
1 5 

2 Founer  components  and  found  the  results  consistent,  signifying  the  sampling 
of  Q(x)  to  be  adequate. 


3.  RESULTS 


’n  this  section  we  compare  the  propagation  effects  of  rod  profile 

striations  and  Gaussian  striations.  Since  the  computer  routine  calculates 
15 

the  signal  at  2 values  of  x over  a 64  km  range,  we  have  a good  sampling 
for  statistics.  The  output  of  the  routine  includes  the  mean  intensity, 
jup,  which  should  be  unity  by  conservation  of  energy  (this  check  on  the 
results  was  unity  for  all  cases),  and  the  scintillation  index, 

S,  -.J\^  . 1 . 

The  distribution  of  the  intensity  is  obtained  by  printing  the  number  of  samples 
in  bins  3 dB  wide.  The  following  table  has  the  percent  of  samples  in  each  bin 
for  those  cases  where  at  least  0.55^  of  the  samples  were  less  than  -3  dB. 

The  columns  begin  with  the  percent  less  than  -30  dB  and  work  up  in  3 dB 
increments  with  the  last  column  for  samples  ^ 12  dB.  The  T signifies  a finite 
percentage  less  than  or  equal  to  .05.  The  calculations  were  carried  out  for 
Ad  up  to  the  lesser  of  10  and  We  chose  to  stop  the  tabulation  at 

Ad  the  lesser  of  10  and 

The  results  for  Gaussian  profiles  do  not  enter  the  table  for 

Ad  < (A*  The  results  for  rod  profiles  do  not  enter  the  table  for 

rms 

Ad  < within  the  range  of  Ad  of  interest. 

If  one  considers  the  question  of  how  one  scales  a critical  quantity 
with  wave  length,  one  will  get  different  scaling  for  the  rod  and  Gaussian 
profiles.  For  example,  consider  the  value  of  on-axis  electron  density,  n^, 
below  which  more  than  99.5^  of  the  signal  is  above  -3  dB.  For  the  Gaussian 
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profile, 


and  n^  scales  as  ,\  for  a given  structure  and  d.  For  the  rod  profiles. 


> 10-V2 


r ^dQ 
e 0 ^rms 


and  scales  as  a~  for  a given  structure  and  d. 

Therefore,  the  important  question  of  the  payoff  of  higher  frequencies 
-2  -3/2 

■:s  bounded  by  the  X and  A ' scaling  for  Gaussian  and  rod  profiles. 

Earlier  wo-'k  by  Mission  Research  Corporation  (MRC)  had  suggested  that 
propagation  effects  depended  on  alone,  which  means  n^  scales  as  . 

Tnis  is  because  their  analysis  assumed  Ad  » 1. 

We  illust.-'ate  this  difference  by  including  in  the  output  the  average 
i;:t  error  probability,  for  noncoherent  frequency  shift  key  (FSK)  modulation 

4 

acc.orriir.'  to  KRC's  formula  , 


= /p:-)?'(T)dY 


where  P(-)  ■:$  toe  probability  density  that  the  signal  to  noise  radio,  (S/N), 

's  T and  PgiY)  -we'  is  the  hit  error  probability  for  a given  (S/N)  of  Y. 

2 

if  Y^^  is  the  (S’N)  in  the  absence  if  the  srriations  then  we  use  Y 3=  '’piUI  and 


/p^(v:iP(Y)dv  = - p; 


where  the  sum  is  over  the  N = 2 values  of  x.  We  illustrate  the  results 
in  Figure  3,  where  the  line  labeled  MRC  is  from  their  report.  The  case 


] 

i 

! 

! 


1 


chosen  is  for  "Y^  = 31.623  {15  dB).  The  bracketed  lines  near  MRC's  curve  show 

2 

our  results  for  Xd  = 10  km  . The  brackets  signify  the  spread  between  rod(r) 

and  Gaussian(g)  profiles.  We  approach  the  MRC  result  at  the  high  limit  of 

Xd.  To  the  right  of  the  MRC  curve  are  three  lines  signifying  the  results 

for  Xd  = O.l(r),  Xd  = 0.1(g)  and  Xd  = 0.01(g).  Sprinkled  in  between  are  various 

other  points  for  other  cases.  It  is  apparent  that  at  lower  values  of  Xd, 

gets  better  faster  for  Gaussian  than  rod  profiles.  The  best  way  to  present 

results  of  this  type  is  on  a Xd  versus  A4  plane.  We  demonstrate  for  two 

rms 

_4  _2 

contours  of  P„  = 10  and  10  in  Figure  4.  The  mesh  points  are  the  cases  we 
have  run.  We  sketch  the  contours  so  they  are  on  the  correct  side  of  the  mesh 
point  values. 

The  utility  of  the  chart  is  apparent  for  putting  the  system- 

2 5 3 5 

effect  problem  in  perspective.  Since  10  £ d £ 10  , in  general  we  see  that 

UHF(x=10  ) lies  in  a band  at  the  top  where  trouble  begins  in  the  region 
0.3  < A$  < 1 for  both  profiles.  S Band(X=10"^)  is  below  this  band  where  trouble 
begins  in  the  region  1 < < 3 for  rods  and  3 < < 20  for  Gaussian 

profiles.  Finally,  K Band(X=10"®)  Is  below  this  band  with  3 < A4>^g  < 10  for 
rods  and  20  < < 100  for  the  Gaussian  profiles.  Thus,  the  highest 

frequencies  are  the  most  affected  by  uncertainties  in  the  profile  form. 


This  chart  was  used  by  researchers  in  an  attempt  to  put  the  measured 
scintillation  index  in  perspective.  We  sketch  contours  of  for  Gaussian 
and  rod  profiles  in  Figures  5 and  6.  The  Gaussian  profiles  give  rise  to  large 

focusing  effects  which  lead  to  large  values  of  in  the  focal  region  mentioned 
in  previous  work^.  From  the  table  we  saw  that  only  Gaussian  profiles  produced 


Signals  >12  dB.  The  rod  profiles  do  not  give  rise  to  high  values  of  S^,  and 
there  does  not  seem  to  be  any  consistency  as  to  which  values  are  above  or 

below  1 as  can  be  seen  by  the  uncertainty  in  drawing  the  = 1 contour  in 
Figure  6. 

We  believe  should  be  abandoned  as  a measure  of  the  severity  of 
propagation  effects  since  it  is  influenced  by  intensity  peaks  which  are 
unimportant  for  degradation.  Comparisons  of  Figures  4.  5 and  6 show  correlation 
Of  the  Ohset  Of  trouble  (p;  = lo’^  and  = 0.3)  but  we  should  have  a better 

measure  of  the  degree  of  trouble.  As  an  exaitple,  we  show  in  Figure  7 the 
contour  for  1%  probability  of  fades  >12  dB  for  rod  and  Gaussian  profiles. 

Here,  the  separate  scaling  is  apparent.  The  correlation  of  this  contour  with 
the  = 0.8  contour  of  Figure  6 is  not  bad  for  a small  range,  but  the  corre- 
lation with  in  Figure  5 for  the  Gaussian  profiles  is  non-existent. 
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ATTN:  William  F.  Ftlaut 
ATTN.  L.  A.  Berrv 

National  Oceanic  .At.mo.-:phene  Admin. 

Environmental  Reseaixh  laboratories 
ATTN:-  Josei)!i  11.  Pope 
ATTN:  C.  L.  Rcienach 

DEPARTME.NT  OF  DEFENSE  CO.NTHACT  Olf^ 

Acromilronic  Ford  Con>oralion 
Western  Development  Lalxiratories  Div. 

ATTN:  J.  T.  Mattinglev,  MS  X22 


UEL’AIU’MLM  or  UM-  KN-SK  CONfHACTUKS  ( Conti iiucj ■ 


! 


i 

I 


Ai  rospiico  CiinH)r;itioi> 

ATTN:  S,  p.  Ik)\ver 
ATTN:  Irung  M.  Garfuiikcl 
ATTN:  SMTA  for  PWW 
ATTN:  T.  M.  Sulmi 
AT  I'N:  V.  Josophaon 

Ana!>tical  Syslcnia  liiigmoeriag  Coiii. 
ATTN:  liasUo  Soicnces 

The  Boeini>  Coiupoiiy 

ATTN:  D.  Murray 
ATTN:  Glen  Keister 

I'mv.  of  California  at  San  Diego 
ATTN:  Henry  G.  Uooker 

C.ilspan  Coiiioration 

ATTN:  HomeoA.  Dcliberis 

Computer  Sciences  Corporation 
ATl'N;  John  Spoor 
ATTN:  11.  Blank 

Comsat  Lalxiratories 

ATl'N:  R.  li.  Taur 


bU’AaT.MkNT  Of  Dkrk.NSK  CUM  RACTOliS  (Continued! 

Institute  for  Defense  Analyses 
ATTN:  Krnest  Bauer 
ATTN:  Hans  WolDtard 
•ATTN:  J.  M.  Aein 
.ATTN:  Joel  Hengston 

Inti.  Tel.  & Telegraph  Corporation 
ATTN:  Tech.  Lib. 

Johns  Hopkins  University 

Aiiplied  Physics  Utlioratory 

ATTN:  Document  Librarian 

Loekliced  Missiles  & ^mce  Co.,  Inc. 

•ATTN:  Dept.  60-12 

Lockheed  Missiles  and  Space  Comjiany 

ATTN:  Billy  M.  McCormac,  Dept.  52-M 
ATTN:  .Martin  Walt,  Dept.  52-10 
ATTN:  Richard  G.  Johnson,  Dept.  52-12 

M.l.T.  Lincoln  Lalxiratorj' 

ATTN:  Mr.  Walden,  Xllii 

ATTN:  Lib.  A-082  for  David  .M.  Towle 

■ATTN:  D.  Clark 

ATTN:  James  H.  Pannell,  L-24G 


Cornell  University 

Department  ot  Klectrieal  Engineering 
ATl'N:  D.  T.  Farley,  Jr. 

ESL,  Ihc. 

ATTN:  James  Marshall 
ATTN:  J.  Roberts 
ATl'N:  V.  L.  Mower 
ATTN:  R,  K.  Stevens 

General  Electrie  Coinpanj 
TE.MPO-Ccntcr  for  Advanced  Studies 
ATTN:  Don  Chandler 
ATTN:  DASIAC 

General  Electric  Companv 
ATTN:  F.  A.  Reibert 

General  Research  Corporation 
ATTN:  John  Isc,  Jr. 

Gcoph)  sical  bistitutc 
Universiw  of  Alaska 

A'l  l'N:  Neal  Brown 
ATTN:  T.  N.  Davis 
ATTN:  Tech.  Lib. 

GTE  Sylvama,  Inc. 

Electronic.s  Sj'stems  Gqi-Eastem  Div. 
ATTN:  Marshal  Cross 

liRB-Singer,  Inc. 

.ATl'N:  Larry  Feathers 

I niver.sitt'  of  Illinois 
Department  of  Electricai  Engineering 
ATTN:  K.  C.  Yeh 


Marlin  Marietta  Corjioration 
Denver  Division 

ATTN:  Special  Projects  Program  2*18 

Maxtvell  Laboratories,  Inc, 

ATTN:  A.  J,  Shannon 
ATTN:  Victor  Fargo 
ATTN:  A.  N.  Rostocker 

McDonnell  Dougla,s  Corporation 
ATTN:  J.  .Moule 
ATTN:  N.  Hams 

Mission  Research  Corporation 
ATTN:  R.  Hendrick 
ATTN:  .M.  Schoibe 
-ATTN:  R.  Bogusch 
ATTN:.  Steven  L.  Gutsche 
ATTN:  P.  Fischer 
ATTN:  Dave  Sowle 

The  .Mitre  Coriioration 

ATTN:  S.  A.  Morin,  .M/S 
ATTN:  C.  E.  Call.-tlnan 
ATTN:  J.  C.  Keenan 
ATTN:  Chief  Scientist,  W.  Sen 
ATTN:  G.  Harding 

The  .Mitre  Corporation 

■ATTN:  .Allen  Schneider 

Pacific-Sierra  Research  Corii. 

■ATTN:  E.  C.  Field,  Jr. 

Photometries,  bic. 

.ATTN:  Ii'Ving  L.  Kofskj 

Physical  Dynamics,  Inc. 

ATTN;.  .a..  Thompson 
-ATTN:.  Joseph  B.  Workman 


R & D Associates 

tvlt'  Gabbard 

ATIN;  Riehard  Latter 

ATTN-  H'h‘Tt“-  "'*‘‘®''‘-  •^‘^■ 
3 cv  A?'m'  ? LeLevIcr 

•«ej  AT  IN:  Forrest  Gilmore 

'ITie  itand  Corporation 
ATTN:  Cullen  Crain 

Science  Applications,  Inc. 

Straker 

ATLN:  D.  Sachs 
ATTN:  Lewis  M.  Unson 

Science  Applications,  Inc. 

Huntsville  Division 

ATTN:  Dale  II.  Dlvis 

•«‘=°rporated 
"t  iiN:  B.  Adams 


Stanford  Research  InstiUite 

Charles  L.  Rino 
ATTN;  G.  Smith 

ATTN 

Nuilson 

At  IN;  Alan  Burns 
ATTN:  Walter  Jayo 

Ai?n  Chesnut 

Aim  L.  L.  Cobb 

ATTN;.  E,  j,  Eremouw 

A-J^n ‘ Corporation 
ATTN.  E.  G.  Meyer 

Tri-Com,  Inc. 

ATTN:  Darrel  Murray 

TRW  Systems  Group 

a;WN:  R.  k.  Plebuch 

AFTN:  Robert  M.  Webb,  MS  Rl-uso 


